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Abstract 


This is the final report for work performed between June 1999 through May 2000 
under contract NAS5-99221. The work represents continuation of the previous contract 
NAS5-32371 which encompasses five areas: 

(1) continued refinements and applications of the 2-D chemistry-transport model (CTM) 
to assess the ozone effects from aircraft operation in the stratosphere; 

(2) studying the mechanisms that determine the evolution of the sulfur species in the 
aircraft plume and how such mechanisms affect the way aircraft sulfur emissions 
should be introduced into global models; 

(3) the development of diagnostics in the AER 3-wave interactive model to assess the 
importance of the dynamics feedback and zonal asymmetry in model prediction of 
ozone response to aircraft operation; 

(4) the development of a chemistry parameterization scheme in support of the global 
modeling initiative (GMI); 

(5) providing assessment results for preparation of national and international reports 
which include the “Aviation and the Global Atmosphere” prepared by the 
Intergovernmental Panel on Climate Change, “Assessment of the effects of high- 
speed aircraft in the stratosphere: 1998” by NASA, and the “Model and 
Measurements Intercomparison II” by NASA. 

Part of the work carried out under NAS-99221 was reported in the final report for 
NAS5-32371. (Note that NAS5-32371 ended in February 1998 and the report was 
written in October 1999. This was done for completeness because many of the reports 
mentioned in (5) were not finalized until that time. This report is written as an update to 
the NAS5-3237 1 report and we will often refer to results reported there.) 

We participated in the SAGE III Ozone Loss and Validation Experiment (SOLVE) 
campaign and we continue with our analyses of the data. 
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I. Ozone Assessment Calculations Using the AER 2-D CTM 

The improvements made to both the chemistry and transport schemes in the AER 
2-D CTM were reported in NAS5-32371. The chemistry scheme now includes 
calculations of effective zonal mean rates for heterogeneous reactions on Polar 
Stratospheric Cloud taking into account temperature variation along longitudes. The 
transport scheme simulates the effect of the tropical barrier through adjustment of Kyy in 
the lower stratosphere. A number of refinements were also made to the 2-D aerosol 
model. The models were used to provide results for the Intergovernmental Panel for 
Climate Change (IPCC) report (IPCC, 1999), NASA’s High Speed Research (HSR) 
Program report (Kawa et al., 1999), and the Model and Measurement Workshop report 
(Park et al., 1999). 

Debra Weisenstein also took the lead to prepare The IPCC Technical Report, entitled 
Model Calculations of the Atmospheric Effects of Supersonic Aviations Which Were Used 
in the IPCC Special Report on Aviation and the Global Atmosphere. The report is 
intended to provide additional detail regarding the supersonic emission scenario 
calculations presented in Chapter 4 of Penner et al. (1999). Due to the space limitations 
in the Special Report, detailed plots and discussion for a majority of the model 
calculations were omitted. The technical report was intended to fill that gap and has been 
compiled for two reasons: (1) to aid those called upon to review the IPCC Special 
Report, and (2) to archive and document the calculations performed by nine modeling 
groups for up to 50 scenarios each. It has been published in electronic form only and 
may be obtained via ftp download in postscript (.ps) or Adobe Acrobat (.pdf) formats at 
the following locations: 

ftp://typhoon.larc.nasa.gov/pub/Sage/IPCC_TECH_REPORTS/supersonic/report 
ftp://ftp.aer.com/IPCC_TECH_REPORT (not for web browsers) 

ftp://anonymous @ ftp.aer.com/IPCC_TECH_REPORT (for web browsers) 

and additionally can be read online via a web browser at: 

http://asd-www.larc.nasa.gov/sage/IPCC_TECH_REPORTS/tech.html 

http://www.aer.com/groups/chem/ipcc_report/. 
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II. Intercomparison of AER, GSFC, and LLNL 2-D Models 

An intercomparison study was performed between the AER, GSFC, and LLNL 2-D 
models in which the transport formulations of those models were run within the AER 
model at their native grid resolutions to simulate the change in ozone in response to a 
fleet of supersonic aircraft. An inert tracer calculation with source similar to HSCT- 
emitted NO y was run to confirm that the transport was implemented properly. Small 
differences were found between the AER model running the GSFC circulation and the 
GSFC model itself due to the different transport schemes used (Lin and Rood for GSFC, 
Smolarkiewicz for AER). It was found that substantial differences between 2-D model 
results still exist even with identical transport, implying that the difference in chemistry 
treatments (particularly PSC treatment) is still significant and should be resolved. 

III. Modeling Activities Analyzing Microphysical Processes in HSCT Wakes 

We developed the AER Far-wake model to simulate the evolution of sulfur species 
in the far-wake and determine how the background aerosol sulfate layer is affected. Our 
results showed that the ozone response to a HSCT fleet is sensitive to the SO 3 /SO 2 ratio 
at the nozzle plane and the possibility of heterogeneous conversion of SO 2 to sulfate in 
the near-field. Changes in sulfate surface area caused by sulfur content in the fuel 
calculated by the AER model were used by all models in the assessment calculations for 
the IPCC report and the HSR report. 

We cooperated with Drs. B. Karcher and R. Busen (DLR Institute of Atmospheric 
Physics, Germany), R.P. Turco and F. Yu (UCLA), and R.C. Miake-Lye (Aerodyne 
Research, Inc.) in estimating global implications of the available particle samplings in 
aircraft wakes. The main results of this cooperation are presented in the Karcher et al. 
(2000) paper, submitted to Journal of Geophysical Research (see Appendix 1). In this 
report, we briefly outline our contribution to this study. 

A phenomenological model has been developed, which allows estimation of the 
aerosol size distribution in aircraft wakes based on such given parameters as fuel sulfur 
content, fraction of emitted sulfur converted to S0 3 , emission index of particulate organic 
matter (POM), wake dilution factor, etc. It was shown that this model could successfully 
explain existing wake sampling data for both subsonic and supersonic aircraft. However, 
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this model works only for a relatively young wake (i.e. less than several hours after 
exhaust). In an old wake, interactions with background aerosol and effects of sulfur 
photochemistry could be important. 

We applied our far-wake (Danilin et al., 1997) and 2-D (Weisenstein et al., 1997) 
models to estimate effects of the different HSCT emission scenarios on the global sulfate 
aerosol layer. Our far-wake model was initialized from the phenomenological model for 
several scenarios and performed calculations up to 2 days. The output of the AER far- 
wake model (aerosol size distribution and sulfur mass partitioning among S0 2 and H 2 S0 4 
in gas and liquid phases) is used to initialize the AER 2-D model. The difference 
between global model runs with and without aircraft sulfur emissions provides the 
response of the global sulfate aerosol layer. The novel approach used in this study 
improves the reliability of the previous global studies (e.g., IPCC, 1999), since it 
constrains them by available wake sampling data. 

IV. AER 3- Wave Model 

The 3-wave model was used to assess the atmospheric effects of HSCT engine 
emissions. The results showed that the feedbacks reduce the ozone depletion in northern 
midlatitudes by 20%, mainly induced by the feedback due to changes in circulation. The 
reader is referred to the NAS5-32371 for details. There is no new work to report here. 

V. Chemistry Parameterization Scheme 

We developed a set of polynomials that will calculate the 24-hour averaged 
production and loss rates of 20 transported species in terms of 18 independent variables. 
The independent variables are temperature, altitude, latitude, time of the year, sulfate 
surface area, PSC surface area, and concentrations of species from the 3-D CTM from the 
previous time-steps. Our test shows that the overall root-mean-square accuracy is about 
10% although differences for individual points may be larger. The details were reported 
in NAS5-32371. This work has been published in JGR in collaboration with Harvard 
University (Duncan et al., 2000). A copy of the paper is included as Appendix 2. 
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VI. Modeling Support for the SOLVE Campaign 

We used funding from this project to perform model analysis of the possible self- 
sampling wake measurements by ER-2 during the SAGE HI Ozone Loss and Validation 
Experiment (SOLVE). The SOLVE campaign provided a unique opportunity for in situ 
measurements of the gas composition and particles in the ER-2 wake under cold 
conditions in the Arctic lower stratosphere during winter. These measurements could be 
crucial for our understanding of the projected HCST fleet impact on the atmosphere, 
since they could provide the necessary constraints for model calculations under PSC-type 
conditions characterized by non-linear chemistry. To justify the importance of this 
experiment and demonstrate the possibility of obtaining valuable measurements with the 
existing ER-2 payload, we performed model calculations for the typical ER-2 exhaust 
scenario and different ambient conditions. Details of our calculations are shown in 
Appendix 3. Indeed, we have shown that the ER-2 payload, despite existing limitations 
of the accuracy and frequency of measurements, could allow us to obtain the necessary 
information about the ER-2 wake composition. 

During the SOLVE campaign, two ER-2 wake interceptions occurred during the 
flights on January 31 and March 7, 2000. The first event on January 31 was 
characterized by weak perturbations of C0 2 and H 2 0, indicating that measurements were 
made at the very edge of the ER-2 wake. However, the second event, on March 7, 
showed a strong increase in NO y , H 2 0, and particle concentrations, signaling that the ER- 
2 wake was penetrated. Our further efforts will be focused on this event. 

According to the SOLVE data exchange protocol, we cannot present these 
measurements in this report. However, we would like to report that we are working 
closely with the Pis of the relevant instruments. Also, we cooperate with Drs. R.C. 
Miake-Lye and R.C. Brown (Aerodyne Research, Inc.), who can provide the necessary 
input data for our far-wake model calculations. Results of our model analysis will be 
presented at the SOLVE Science Team Meeting in September 2000. 
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VII. Various Meetings 

VII.l Meeting on Next Generation Supersonic Civil Transport Risk 
Assessment Methods, Washington DC, November 1999 

This is a workshop organized by EPA and NASA to advise ICF Consulting on 
their plan to assess the health risks associated with oozne depletion that may arise from 
aircraft operation. The meeting took place at the Office of ICF Consulting, Washington, 
DC, November 18-19, 1999. Malcolm Ko was invited to the meeting and participated in 
the preparation of the workshop summary. A second meeting is being scheduled in 
August 2000. 

Vn.2 GMI Meeting, Washington DC, January 2000 

Malcolm Ko attended the GMI meeting January 2000 and made a presentation on 
‘The possible future roles for 2-D models in GMI”. The presentation was made in 
collaboration with Charles Jackman from GSFC. A copy of the viewgraphs is included 
in Appendix 4. 

VII.3 AESA Meeting, Snowmass CO, June 2000 

Debra Weisenstein and Michael Danilin attended the AEAP workshop at 
Snowmass CO. A GMI meeting was scheduled. Michael Danilin made a presentation 
entitled “Global Implications of the Ozone Loss in a Space Shuttle Wake”. A copy of the 
viewgraphs is included in Appendix 5. 
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On the Unification of Aircraft Ultrafine Particle Emission Data 


B. Karcher, R.P. Turco, F. Yu, M.Y. Danilin, 
D.K. Weisenstein, R.C. Miake-Lye and R. Busen 


A-l 
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On the Unification of Aircraft Ultrafine Particle Emission Data 


B. Karcher 1 , R.P. Turco 2 , F. Yu 2 , M.Y. Danilin 3 , D.K. Weisenstein 3 , R.C. Miake-Lye 4 , R. Busen 1 


Abstract 

To predict the environmental impacts of future commercial aviation, intensive studies have been 
launched to measure the properties and effects of aircraft emissions. These observations have 
revealed an extremely wide variance with respect to the number and sizes of the particles 
produced in the exhaust plumes. Aircraft aerosol ultimately contributes to the population of 
cloud-forming nuclei, and may lead to significant global radiative and chemical perturbations. In 
this paper, recent discoveries are coordinated and unified in the form of a physically consistent 
plume aerosol model that explains most of the observational variance. Using this new approach, 
it is now practical to carry out reliable global atmospheric simulations of aircraft effects, as 
demonstrated by a novel assessment of the perturbation of the stratospheric aerosol layer by a 
supersonic aircraft fleet. 
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1. Introduction 

Aviation impacts atmospheric aerosols and clouds 
by emitting copious amounts of small liquid aerosol 
and soot particles [ Seinfeld , 1998; Boucher , 1999; In- 
ternational Panel on Climate Change ( IPCC ), 1999]. 
However, a quantitative assessment of potential ef- 
fects requires a better understanding of ultrafine (size 
< 10 nm) particle formation mechanisms. Besides 
sulfuric acid (H2SO4), charged molecular clusters 
(chemi-ions) and condensable organic species have 
been suggested to explain important aspects of the 
aerosols observed in jet aircraft plumes [ Yu and Turco , 
1997; K archer et al 1998a]. All of these species have 
now been measured in-situ under cruise conditions 
[<7iirttU5 et al ., 1998; Slemr et al., 1998; Arnold et al., 
1999]. 

To integrate such field data with physical theory to 
create a more useful model for aircraft emissions, a pa- 
rameterization scheme is derived here that accurately 
predicts the properties of volatile jet aerosols. The 
model is consistent with microphysical mechanisms 
describing the phenomenology seen in aircraft plumes, 
including processes driven by electrical charge and or- 
ganic vapors. We show that existing plume particle 
measurements obtained over a wide range of flight 
conditions can be unified within the framework of our 
new parameterization, dramatically reducing the ap- 
parent variability in volatile particle emission indices 
previously inferred from the field data. 

The emission model is finally applied to carry out a 
new globed assessment of aircraft impacts on ambient 
aerosol properties, narrowing considerably the range 
of projected enhancements in background aerosol sur- 
face area densities associated with commercial air 
traffic. Hence, we can demonstrate that this approach 
can be employed in future evaluations of atmospheric 
perturbations caused by aviation. 

2. Lessons From Field Observations 

In addition to emitting nonvolatile soot particles 
formed within the engine combustors, subsonic and 
supersonic jet aircraft generate volatile particles within 
their near- field exhaust plumes. Models [e.g., Danilin 
et al., 1998] and observations [e.g., Hofmann et al. , 
1998] show that the total background aerosol num- 
ber and surface area densities can be significantly en- 
hanced by aircraft emissions. The resulting tenden- 
cies of aircraft sulfate and soot aerosols to influence 
cirrus clouds and to alter chemical processes are as 
yet poorly quantified, limiting our ability to assess 


climate and ozone impacts [Kdrcher, 1999]. Indeed, 
in-situ observations of particles in fresh (plume age 
t < Ihr) aircraft plumes show a considerable scat- 
ter in the apparent volatile particle emission index, 
even for a fixed fuel sulfur (S) emission index EI(S) 
[IPCC, 1999; Anderson et al., 1998; Schroder et al., 
1998; Brock et al, 2000]. The factors that contribute 
to this variability have not been indentified. 

It had been proposed that new volatile aerosols 
are composed mainly of aqueous H2SO4 [Hofmann 
and Rosen, 1978]. However, some measurements 
do not support this presumption [Miahe-Lye et al, 
1998; Kdrcher et al, 1998a; Yu et al, 1999; Brock 
et al, 2000]. Simulation models that include con- 
current emissions of condensable organic species re- 
sulting from incomplete fuel combustion along with 
chemi-ions also formed in the combustors can ex- 
plain recent field observations for various S emission 
levels [Kdrcher et al, 1998a; Yu et al, 1999]. In 
contrast, models that neglect organic emissions and 
chemi-ions either fail to explain observations or re- 
quire H2SO4 emissions exceeding values measured in- 
situ [Curtins et al, 1998] or calculated using chemi- 
cal kinetic codes [Lukachko et al, 1998; Tremmel and 
Schumann, 1999]. Nevertheless, the emission indices 
of organic compounds are not yet well established. 

3. Unified Representation of Volatile 
Particle Properties 

3.1. Model Description 

We have developed a conceptual model to explain 
observed ultrafine particle emissions that consolidates 
several new concepts concerning the physico-chemical 
mechanisms responsible for the nucleation and growth 
of aviation plume aerosols. With this model, we 
have carried out sensitivity studies to identify the 
key parameters governing particle formation. Con- 
sistent with field data, we assume that particulate 
organic matter (POM) and H2SO4 initially condense 
to form both ionized (subscript z) and electrically 
neutral (subscript n) aerosol modes. This phase of 
particle formation occurs at plume ages shorter than 
0.1s. The effects of electrical charge on the particle 
collection kernels lead to a preferential growth of an 
‘ion’ aerosol mode, whereas the neutral mode particles 
maintain smaller cluster sizes and are not detectable 
by condensation nuclei counters (CNCs). Scavenging 
of the neutral clusters by ion mode particles consti- 
tutes the major growth mechanism of the latter dur- 
ing the first minutes of plume evolution. Restrict- 
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ing the analysis to such times allows the interaction 
of volatile particles with emitted soot and entrained 
ambient aerosol to be neglected. 

The phenomena noted above can be expressed by 
means of two simplified equations: 

EI(POM) + ^EI(S) = f E PkWkN k D\ (1) 

k=i>n 

^ = -KmN n . (2) 

Equation (1) describes the overall mass conservation 
between the prescribed emissions of POM and H2SO4 
(in units of g emitted per kg of fuel burned, where 
77 is the conversion fraction of fuel S (32 g mol" 1 ) to 
H2SO4 (98 g mol" 1 ) molecules at emission) and the 
particle mass contained in the two volatile modes (p 
and W are the mass density and mass fraction, re- 
spectively, of POM or H2SO4). Assuming log-normal 
size distributions, N and D denote the total particle 
number per kg-fuel and the volume mean diameter, 
respectively, for each mode. Beyond a timescale of 
about 0.1s, Ni is constant and equal to the number 
of chemi-ion-induced aerosol particles, which is con- 
trolled by ion-ion recombination at earlier times [Yu 
and Turco , 1997]. Equation (2) describes the rate of 
change in the number of neutral clusters, where the 
coagulation kernel K depends on D,. A kernel for 
thermal coagulation in the free molecular regime of 
the form K = K 0 (Di/D i0 ) 2 is used, where D i0 = 
1 — 4nm is the initial value (depending on EI(POM) 
and on 7 ] EI(S)) and Kq = 1.2 x 10" 9 cm 3 s _1 is the 
appropriate kernel at that plume age. The choice 
K ( x D 2 approximates the collection kernel of small 
neutral clusters by the growing ion mode particles. 
(Turco and Yu [1999] present the analytical solutions 
of surface- controlled collection of particles in an ex- 
panding plume, among other cases, in a more general 
context.) D n is assumed to be constant and equal to 
1 nm, since the growth rate of the neutral clusters is 
slower than of the ionized particles. The plume di- 
lution history enters the model through n^, the num- 
ber density of ion mode particles per cm 3 of air. In 
(2), rij = CNi'D(t) 1 where the plume dilution factor 
[Schumann et a/., 1998] is taken to be V = (t m /t) a 
(£ m = 1 s and a — 0.8) and where C ~ 6 x 10" 11 kg- 
fuel cm" 3 (evaluated at subsonic cruise conditions) 
depends on the mass density of air. 

Equations (1) and (2) can be solved analytically to 
obtain the temporal evolution of Di(t). We simplify 
aerosol thermodynamics and set p — 1.5g/cm 3 and 
W = 0.6 for both modes. The maximum ion mode 


volume mean diameter 

_ r EI(POM) + EI(H a S0 4 ) l */ 3 
Di ~ l npWNi/6 J ( ' 

is approached within plume ages t ~ 10 3 -10 4 s and lies 
within the range 6-12 nm. Equation (3) follows di- 
rectly from (1) with EI(H 2 S0 4 ) = 98t ? EI(S)/ 32. Dif- 
ferentiating (1) and then substituting from (2) allows 
the temporal evolution of Di(t) to be derived in an- 
alytical form. In terms of the normalized variables 
x = Di/Df° and r = t/t m we find: 



Here k = KoNiCt m (Di° / Dio) 2 /3 is the single dimen- 
sionless parameter governing the solution obtained in 
the form r(x), from which the temporal evolution 
of Di(t) can be derived. The left-hand integral is 
solved by the method of partial fractions, yielding 
F(x) - F(x 0 ) y with F(x) = (1/6) ln[(l + x + x 2 )/(l - 
x) 2 ] -1- (1/v^) arctan[(l + 2x)/\/3]. The right-hand 
integral yields (r 1 "" — TQ -ot )/(l — a) for a / 1 or 
ln(r/ro) for a = 1. 

To facilitate comparison with observations, we de- 
fine the apparent particle emission indices, AEIs, as 
the cumulative number of particles per kg-fuel that 
are larger than a given size. Hence, the AEI is the con- 
volution of the size distribution of the ion mode par- 
ticles inferred from (4), after the detection efficiency 
of the CNC has been taken into account. Although 
great effort has been spent in calibrating CNCs in 
the laboratory [e.g., Cofer et al , 1997; Schroder et 
al, 1998; Brock et al , 2000], their performance under 
flight conditions remains poorly defined. To incorpo- 
rate the effects of variable CNC detection efficiencies, 
e, to a first order and to permit an analytical evalua- 
tion, we assume a linear response function 

e(D,D th ) = 0.5 + ^-=^, (5) 

where 5 is the width of e, that is, e = 0 for D < 
Dth — S/2 and e = 1 for D > D t h + S/2. We use 
S ~ D t h as a free parameter. To obtain the ap- 
parent emission indices AEI, we integrate e over the 
log-normal size spectrum F(D, y a), correspond- 
ing to the ion mode particles, assuming a modal 
width a = 1.3 and using the mean number diam- 
eter Df*(t) = D{(t) exp(-31n 2 a) deduced from our 
analytic solution: 

AEI(t) - Ni r dD F(D, D?,o)e(D, D th ) . (6) 
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The choice of log-normal size distributions with a = 
1.3 is supported by a detailed comparison of numer- 
ical simulations with near-field data [Karcher et al, 
1998a; Yu et al , 1999], 

3.2. Discussion of Model Features 

Figure 1 illustrates the behavior of the calculated 
AEIs as a function of EI(H2S04) oc r;EI(S), for vari- 
ous plume ages (la), POM Els (lb), and lower CNC 
detection limits (lc). Common to all of the results in 
Figure 1 is the decreasing sensitivity of AEI for H2SO4 
levels above 50— 100 mg/kg fuel (or ~ 0. 3-0.6 g S/kg- 
fuel using 77 = 0.05; the global average S content is in 
the range 0.4— 0.6 g S/kg-fuel). In Figure la, the rapid 
increases of AEIs at smaller H2SO4 levels are caused 
by the growth of ultrafine particles across the detec- 
tion threshold D t h over the course of time. In con- 
trast, for high H2SO4 levels, growth takes place at ear- 
lier stages (£ < 1 s), so that almost all of the ion mode 
particles (2 x 10 1 7 /kg- fuel in this case) are detected at 
the plume ages indicated. Since AEIs at even lower 
EI(H 2 SC>4) are determined by the condensation of or- 
ganics, the AEIs are very sensitive to EI(POM) in 
this parameter region (Figure lb). The impact of 
CNC cut-off sizes is illustrated in Figure lc. In the 
near-field plume, the use of instruments with rela- 
tively large cut-off diameters (D t h > 5 nm), like those 
employed in the earliest plume measurements, virtu- 
ally precludes the detection of new plume aerosols, 
although the strong sensitivity of AEI on D t h de- 
creases with increasing EI(H 2 SC>4) (Figure lc) and 
plume age (not shown). The behavior of AEIs indi- 
cated in Figure 1 is supported by detailed numerical 
simulations of plume microphysical processes [ Yu and 
Turco, 1997; Karcher et al, 1998b; Yu et al , 1999], 
which increases our confidence in the present repre- 
sentation. 

In comparing the predictions of (1) and (2) with 
CNC data from recent in-situ measurements [Ander- 
son et al , 1998, 1999; Schroder et al, 1998, 2000; 
Petzold et al, 1999], it is noteworthy that the key pa- 
rameters have quite different ranges of variability. For 
example, POM emissions are not well characterized 
observationally while 77 exhibits a range of 0.005—0.15 
andiVj = (1— 4)x 10 17 /kg-fuel. Thus, EI(POM) intro- 
duces the largest uncertainty into our model, but only 
at low H 2 S04 emissions. 

3.3. Comparison With Field Observations 

Figure 2 depicts the dependence on plume age of 
observed (symbols) and modeled (curves) AEIs, for 


the ATTAS aircraft in Figure 2a and for the F-16 air- 
craft in Figure 2b. While the subsonic ATTAS uses 
old- technology jet engines, the F-16 is a military air- 
plane capable of supersonic flight and equipped with 
low-bypass turbofan engines. The F-16 engine com- 
bustor and turbine are similar to commercial high- 
bypass engines and should be representative of re- 
cent commercial technology. In Figure 2a, each sym- 
bol actually represents an average over many individ- 
ual observations, with am associated spread in AEI 
of 30-50% [Schroder et al, 1998]. In Figure 2b, a 
subset of data points [Anderson et al, 1999] is used, 
including the non-volatile soot emissions. All of the 
data points include variations related to differences 
in engine power settings, the location of the measure- 
ment within the plume cross section, and ambient 
conditions (pressures of 220— 360 hPa, temperatures 
of 220— 245 K, and relative humidities of 20—80%). 
Figure 2 demonstrates the ability of the phenomeno- 
logical model to explain the near-field observations 
within the uncertainties of the measurements through 
the selection of reasonable values for EI(POM), 77, and 
Ni. 

The model reflects the generally observed increase 
in AEIs with rising EI(S). A central conclusion from 
this work is that much of the scatter in AEIs noted 
in field measurements can be explained consistently 
by variations of t , EI(POM), and Dth • The seemingly 
pronounced sensitivities noted earlier are due mainly 
to the steepness of the large-size tail of the ultrafine 
particle spectrum. The sensitivity of AEIs to other 
parameters such as the plume dilution rate and am- 
bient conditions is much less pronounced. This does 
not hold when contrails form, because water and ice 
condensation dramatically change the size distribu- 
tions of the plume particles. Due to the large avail- 
able surface area of contrail ice crystals, the number 
of ultrafine particles decreases by coagulation scav- 
enging, even in short-lived contrails [Anderson et al, 
1998; Schroder et al, 1998; Karcher et al, 1998b]. 
These contrail processing effects are not included in 
the present model. 

4. Unifying Plume Particle 
Measurements 

The data discussed above are replotted in Figure 3a 
together with measurements obtained in other field 
missions. Multiple data points for a specific aircraft 
at a given EI(S) correspond to different plume ages 
and CNC cut-off sizes. To arrive at a more unified 
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data set, each point was fitted to the model discussed 
in Section 3.1, obtaining curves similar to those shown 
in Figure 2. Then, each AEI was recalculated using a 
common cut-off size (5nm) and POM emission index 
(20 mg/kg-fuel) at a specific plume age (3 s). Accord- 
ing to Figure 1 , these parameters lead to the largest 
sensitivity in AEI. The resulting normalized data set 
is shown in Figure 3b (symbols). We overlay solu- 
tion curves corresponding to the same £, D t h > and 
EI(POM) at fixed values of 77 and JV*. These compar- 
isons delineate the common trends that underly the 
various measurements made behind different aircraft, 
at different plume ages, and under quite different ex- 
perimental conditions. Given the spread between the 
individual data points (vertical bars in Figure 3a), 
Figure 3b demonstrates that a very wide range of 
observations can be quantitatively explained using a 
rather simple conceptual model with only three basic 
parameters. 

The principle trend of the normalized AEIs with 
EI(S) can be understood even with a constant S con- 
version fraction 77 within the range 0.005-0.05, im- 
plying that 77 varies much less than previously con- 
jectured between different types of aircraft [ IPCC \ 
1999]. However, for the low-S cases, uncertainties in 
the choice of EI(POM) do not allow us to draw a 
definitive conclusion on this point. 

We note two exceptions. Assuming 77 = 0.15 
as perhaps an upper bound, we require a compara- 
tively high value of 70 mg POM/kg-fuel to explain 
the Concorde data, which were performed in much 
more aged plumes (15 min— 1 hr), exhibiting a much 
larger spread in AEI [Fahey et ai, 1995]. While these 
values still fall within the ranges reported in the liter- 
ature [Climatic Impact Assessment Program (CIAP), 
1973; Fahey et ai, 1995; Lukachko et ai , 1998; Trem- 
mel and Schumann , 1999], this case deserves specific 
attention in the future. To explain the F-16 data, we 
assumed negligible organic emissions and low S con- 
version (see Figure 2b). Whether this indicates the 
potential for a substantial variability in the organic 
emissions, or merely reflects a specific feature of this 
type of engine and the fuel used, remains open at 
this time. It should be remarked that discrepancies 
in the determination of EI(S) for the F-16 have been 
reported [Anderson et ai , 1999]: the cause for the dis- 
crepancy between the value of 0.146 g S/kg-fuel used 
here and 0.018 g S/kg-fuel resulting from two inde- 
pendent laboratory analyses is unresolved. 


5. Global Impact of Aviation 
Particulates on Stratospheric Aerosols 

The evaluation of aviation- produced aerosol im- 
pacts on the global atmosphere requires the use of 
a global-scale tracer model in conjunction with a suf- 
ficiently detailed characterization of the aircraft emis- 
sions. To achieve this goal, we have coupled our an- 
alytical formulation for the near-field plume to a far- 
field plume model [Danilin et ai, 1997] and a two- 
dimensional ( 2 -D) global aerosol model [ Weisenstein 
et ai , 1997]. 

The approach used here combines aerosol forma- 
tion in the near-field plume, aerosol processing in the 
far-field plume, and the resulting fleet aerosol source 
and interaction of aerosol processes with atmospheric 
dynamics on the global scale. Our parameterization 
scheme provides unified aerosol size distributions pre- 
scribed according to EI(S), El (POM), and 77 at plume 
ages 10 3 — 10 4 s, excluding the more extreme F-16 
and Concorde cases. The far-field model continues the 
chemical and microphysical processing of the exhaust 
products within the region of greatest aircraft fuel 
burn in the Northern Hemisphere (pressure of 65 hPa, 
temperature of 220 K) for two days as the plume ex- 
pands. The partitioning of S between aerosols and gas 
phase species and the size distribution generated af- 
ter two days are used along with the projected aircraft 
fuel usage for the year 2015 [IPCC, 1999] to obtain 
the aircraft emission source function for the 2-D sim- 
ulations, which assume surface sources of S only for 
the calculation of the background aerosol layer. The 
far field and global models employ the same aerosol 
microphysics to track the aerosol evolution. In those 
models, POM is not dealt with explicitly; instead the 
mass of emitted POM is incorporated into the aerosols 
assuming the same properties as H 2 SO 4 . This as- 
sumption is reasonable, because POM is likely to be 
very soluble in aqueous H 2 SO 4 and may also react 
with the acid to form stable compounds [Kdrcher et 
ai , 1998; Yu et ai , 1999]. 

Figure 4 gives the predicted changes of the am- 
bient sulfate aerosol surface area density (SAD) due 
to the operation of a supersonic aircraft fleet for the 
region of maximum perturbation as a function of 77 
at emission. The simulations assume 0.4 g S/kg-fuel 
and either 0, 20, or 60 mg POM/kg-fuel. The re- 
sults show considerable sensitivity to 77 , with ASAD 
varying from 0.5— 0 . 871 m 2 cm' 3 . The changes in SAD 
should be compared with the simulated background 
SAD of 0.9 jmi 2 cm" 3 . In the stratosphere, ASAD in- 
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creases with rising r] because the plume particles be- 
come larger and are therefore less efficiently scavenged 
by background particles. Varying r\ from 0.005 to 0.1 
implies values of 6 to 120 mg H2S04/kg-fuel, which 
explains the low sensitivity of the results to EI(POM) 
and its decreasing impact when r\ rises. Reductions 
of EI(S) or increases of El (POM) would lead to a 
greater influence of the POM emissions. The present 
simulations are compared with those adopted by the 
IPCC (diamond in Figure 4), which were based on 
the same 2-D model but neglected plume processing 
of the aerosols. The resulting ASAD is close to, but 
slightly larger than, the values plotted for the current 
case, implying that the IPCC calculations are fairly 
consistent with the results derived here. The IPCC 
supersonic scenarios also used higher values of rj (0.5, 
1) for a sensitivity analysis inasmuch as earlier studies 
had not precluded such a range. 

The present analysis indicates that the most im- 
portant parameters affecting large-scale cloud and 
chemical impacts are tj and EI(POM). The reduction 
of both of these factors reduces the particle size and 
thus their atmospheric residence time. Accordingly, 
as technological measures are taken to reduce the sul- 
fur in aviation fuel, the importance of POM will in- 
crease, leading to the need for better information on 
POM emissions, speciation, and heterogeneous chem- 
ical reactions. 

6. Summary and Conclusion 

We have developed a unique unified representation 
of the properties of volatile particles emitted by jet 
aircraft. Our new parameterization is based upon de- 
tailed microphysical mechanisms that incorporate the 
latest phenomenology observed in aircraft plumes, in- 
cluding the roles of chemi-ions and organic vapors. 
The main findings of this study are: 

1. The scatter in emission indices of volatile particles 
seen in the field measurements is due mainly to vari- 
ations of plume age, condensable organic emissions 
(when low S fuel has been used), and the detection 
threshold size of the ultrafine particle counters. 

2. The principle trend of the volatile particle concen- 
trations with fuel S content can be explained with 
conversion fractions of S to H2SO4 at emission within 
the range 0.005 to 0.05, implying that S conversion 
varies much less than previously thought. 

3. The key parameters affecting the large-scale cloud 
and chemical impact of aviation particulates are the 
S-to-H 2 S0 4 conversion fraction and condensable or- 


ganic emissions, the latter especially below H2SO4 
emissions of 50—100 mg/kg- fuel. 

Here, the model has been extensively compared 
with a diverse set of measurements taken in jet plumes, 
and shown to reproduce quite accurately the aerosol 
properties crucial to environmental assessments. Re- 
markably, the complex evolution of the particulates 
can be forecast using only three parameters associated 
with aircraft engines: the amount of emitted chemi- 
ions, sulfuric acid, and condensable organic species. 
We have given uncertainty ranges of these parame- 
ters, within which existing field data can consistently 
be explained with our approach. While chemi-ion and 
sulfuric acid emissions are constrained by field data, 
the a priori unknown speciation and emission levels 
of organic species introduce the largest uncertainties, 
but only at low H2SO4 emissions. 

We have also demonstrated, within the context of a 
global tracer model, that our unified approach is use- 
ful in carrying out detailed simulations of the atmo- 
spheric impacts of aircraft fleet operations. Accord- 
ingly, we have updated and validated the calculations 
used in the preliminary IPCC assessment of aircraft 
effects on climate, and set the stage for future research 
activities. 
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Figure 1 . Apparent ultrafine volatile particle emission indices versus H2SO4 emission index (lower abscissa) and S 
emission index (upper ascissa). Default values used in the analytical model (solid curves) are: iV» = 2x 10 17 /kg-fuel, 
r] = 0.05, D t h = 3nm, EI(POM) = 20 mg/kg-fuel, and t = 3s. a: Plume age t is varied from 0.5 s up to lhr. b: 
Emission index of condensable organic species EI(POM) is varied from 0-80 mg/kg-fuel. c: Lower cut-off diameter 
Dth (50 % detection efficiency) of ultrafine particle counter is varied from 2 — 7 nm. 
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Plume age t, s 

Figure 2. Observed ultrafine particle emission indices (symbols) versus plume age. The legends denote CNC 
detection limits and fuel S contents. The corresponding curves are results from the model described here, a: 
Data taken during the DLR SULFUR 5 and 6 field missions behind the ATTAS aircraft. For the calculations, 
we use 7] ~ 0.03-0.06 depending on EI(S) taken from a detailed chemical model [Brown et al , 1996] and JV* = 
(1.5—2) xl0 17 /kg-fuel [Schroder et al . , 1998]. We used the basically unknown parameter EI(POM) within the range 
20— 30mg/kg-fuel to fit the data. Only volatile AEIs are shown, b: Same as in Figure 2a, but for the F-16 aircraft. 
Data taken during the NASA Subsonic Assessment Near-Field Interaction Flight (SNIF III) Experiment [Anderson 
et al, 1999]. Model curves obtained with r) = 0.005 — 0.01, Ni = (3—4) x 10 17 /kg-fuel, and EI(POM) = lmg/kg- 
fuel. The SNIF data include non-volatile soot particles which show comparatively little variability. Therefore, the 
soot AEIs added to the volatile AEIs create a baseline value of 2 x 10 15 particles per kg-fuel. 
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Figure 3. Emission indices of detectable volatile particles versus fuel S content. Aircraft were examined as part of 
several missions in Europe and the United States of America. Aircraft types and values for D t h are indicated in the 
legend. Besides the data shown in Figure 2, we include a subset of data (those where no contrails formed during 
the observations) corresponding to plume ages up to 2 min (Concorde: up to lhr) discussed by the IPCC [1999]. 
a: Data plotted irrespective of plume ages, CNC detection limits, and organic emissions. The vertical bar at the 
rightmost data point indicates the approximate spread of most of the subsonic data. For the Concorde (uppermost 
triangle and bar), this variability is much larger [Fahey et a/., 1995]. b: Same data as in Figure 3a, but normalized 
to t = 3 s, D t h = 5nm, and 20 mg POM/kg-fuel, after being fitted to (1) and (2) as in Figure 2. The curves 
emphasizing the trend of observed AEIs with increasing EI(S) are model results at the same values for i, D t h > and 
EI(POM), with N{ — 2 x 10 17 /kg-fuel and 77 = 0.05 (solid curve) or 17 = 0.005 (dashed curve). Varying EI(POM) 
and Ni causes the curves to shift, as indicated by the arrows; the resulting shifts are not strictly proportional to 
EI(S). 
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S conversion fraction, r\ 

Figure 4. Calculated changes in sulfate aerosol surface area density in the region of maximum perturbation 
(32-90°N, 14- 18.5 km altitude) versus S conversion fraction at the engine exit for a future (2015 scenario) super- 
sonic fleet. The perturbation should be compared with the simulated background value of 0.9 ^m 2 cm -3 . Curves 
represent calculations with our novel approach, using 0, 20, and 60 mg POM/kg-fuel. The diamond represents 
the simplified treatment used in the IPCC [1999] assessment, repeated for 0.4 g S/kg-fuel and assuming yj = 0.1 
and monodisperse particles with D = 20 nm as the input for the global model, in lieu of aerosol processing on the 
scale of the plume. The diamond is plotted by estimating the engine exit plane value of rj which would give tj = 
0.1 after plume processing. Gas phase chemical conversion of emitted SO2 to H2SO4 and subsequent condensation 
adds about 0.02 per day to 17. 
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Abstract We present a parameterization for the tropospheric concentration of the hydroxyl 
radical (OH) which can be used to overcome the costs of solving kinetic equations in chemical 
tracer models. This parameterization accurately represents OH predicted by a full chemical 
mechanism. The 24-hour average concentration of OH is represented as a set of high-order 
polynomials in variables such as temperature, latitude, declination, and the concentrations of 
ozone, water vapor, carbon monoxide, nitrogen oxides (as a family), and hydrocarbons. The 
parameterization of OH consists of computer- written FORTRAN functions for an efficient 
computation of the polynomials. The parameterization of OH is publicly available. 


1. Introduction 

Three-dimensional (3-D) models of dynamics and chemistry 
have become an important tool for investigating the budgets of 
trace species such as ozone (O3), nitrogen oxides (NO, = NO + 
N0 2 + 2N 2 0 5 +N0 3 + HN0 2 + HNO4), hydroxyl radical (OH), 
and methane (CH 4 ) [e.g.. World Meteorological Organization 
(WMO), 1999]. As computing power has increased, these models 
have become more complex, particularly in terms of the detail of 
the chemical schemes, and have moved toward higher vertical and 
spatial resolution. Many initial studies relied on short (1-2 year) 
simulations, but several problems of interest require longer simu- 
lations, for example, examination of the causes of trends in trace 
gases. Despite the increases of computing power, parameteriza- 
tions of model processes designed to enhance computational effi- 
ciency are still necessary. One such parameterization, that of the 
chemistry of the OH radical, was developed in the late 1980s by 
Spivakovsky et al. [ 1 990a). 

When the primary sink of a species is reaction with OH, a 
parameterization of OH can be used to overcome the expense of 
solving chemical kinetic equations. Spivakovsky et al. [1990a] 
presented a method used to define a set of polynomials describing 
the functional dependence of tropospheric OH on chemical, radia- 
tive, and meteorological variables. Spivakovsky et al [1990b] 
then applied this method to a 3-D simulation of methylch loro- 
form. Using the same method, Jacob et al. [1993] and Wang et al 
[1998] parameterized photolysis rates and chemical production 
and loss rates of various tracers such as O3, OH, and NO,. 
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The parameterization of OH by Spivakovsky et al [1990a] is 
out of date because of significant changes recommended recently 
for key reactions affecting OH concentrations, including a non- 
negligible quantum yield for CK 1 /)) at wavelengths between 312 
to 320 nm [Talukdar et al, 1998) and a decrease in the rate of 
reaction of OH with CH 4 [Gierczak et al, 1997]. In this paper we 
present an updated version of the parameterization of OH pro- 
duced with a current and detailed mechanism of O3, NO, and 
hydrocarbon chemistry [Jaegle et al, 1999]. We include the 
dependence of OH on isoprene and several other hydrocarbons 
omitted by Spivakovsky et al. We improved the method of Spiva- 
kovsky et al. by further automating the parameterization proce- 
dure. The parameterization of OH developed here is intended for 
use in studies of factors controlling distributions and temporal 
trends of tropospheric carbon monoxide (CO) and CH 4 , but can 
also be applied to studies of other chemical species where reaction 
with OH is the primary sink. 

2. Generation of the Parameterization 

Spivakovsky et al [1990a] created their polynomial functions 
by calculating the least squares fit (LSF) to results of a large num- 
ber of point-model calculations. A 0-D photochemical model 
with an accurate kinetic solver was used to generate the point- 
model calculations. These calculations were performed across the 
entire multidimensional domain defined by ranges of independent 
variables, such as 0 3 , CO, water vapor, and NO,. Computer-gen- 
erated FORTRAN functions for evaluation of the polynomials 
were returned as result of this procedure, as well as estimates of 
the accuracy achieved for an independent set of points. If the 
accuracy was found to be insufficient, the domain space was 
divided manually (i.e., off-line) into rectangular subdomains, and 
the parameterization procedure was repeated. (The domain space 
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is divided because at a certain stage in the selection process, aug- 
mentation of additional terms no longer improves the accuracy 
significantly.) Divisions of the domain into subdomains were cho- 
sen based on the knowledge at that time of the dependence of OH 
on various independent variables, such as water vapor, O 3 , and 
NO f . The choice of divisions within a subdomain was often an 
educated trial-and-error process because it is nearly impossible for 
a human mind to conceptualize the dependence of a function on a 
large number of independent variables. 

In an effort to make the parameterization procedure less labor- 
intensive and more efficient, we implemented a computerized 
search for optimal divisions within a subdomain, proposed (but 
not implemented) by Spivakovsky et ai [1990a]. This automated 
procedure is more efficient because it searches for divisions across 
the domain space (i.e., multidimensional divisions). While divi- 
sions with respect to only one independent variable will remove 
some of the nonlinearity in a system, the best divisions for OH 
chemistry usually lie across the subdomain space and involve cor- 
relations of independent variables [Sillman et al. , 1990]. As a 
result, fewer divisions are required to achieve a desired accuracy. 
The division procedure is outlined below. 

The N-dimensional domain is divided into two parts by a plane 

N 

£ a i x i + a N+\ = 0 (1) 

1 = 1 

where denote parameters. Coefficients a, (for i = 1, ... , jV+ 1) 
are obtained as a solution to the following optimization problem. 
Define function p(aj, ... , a/v+i) as a root-mean- square error 
(RMS) for a piece-polynomial function composed of two second- 
order polynomials with a separate LSF polynomial for each of the 
two parts of the domain: 


Table 1. Independent Variable Ranges Used in the 
Parameterization of OH 


Independent Variables 

Ranges 

Pressure, mbar 3 

800-1020, 350-800, 100-350 

Temperature, K 

depends on season, altitude, latitude 

Chemical Variables 


NO„ pptv a 

1-300, 300-1000, >1000 

O 3 , ppbv 

3-110 

CO, ppbv 

10-400 

CH 4 , ppbv b 

600-2200 

H 2 0, ppmv 

depends on season, altitude, latitude 

Isoprene, pptv* 

0-5000 

Acetone, propane, ethane. 

depends on season, altitude, latitude 

propene, >C 4 alkanes, pptv 


Solar Irradiance Variables 


Surface Albedo (unitless) 

0.05-0.8 

Declination Angle, deg a 

'23-'8/8-22, 8-23/21-8 

Latitude, deg a,c 

0-30, 30-40, 40-60, 60-90 North and 
South 

Cloud Albedo (unitless) d 

0 - 0.6 

Overhead O 3 Column, DU 

200-400 

DU, Dobson units. 
a Manual divisions in range (i.e. 

, done prior to parameterization). 


b The range for CH 4 is expanded beyond what is expected in the current 
atmosphere because the present parameterization of OH is intended for 
use in simulations of its past, present, and future tropospheric distribu- 
tions. 

divisions based on resolution of 0 3 climatologies of Logan [1999a,b], 
d Cloud albedo has above and below components. 


Table 2. RMS Errors Associated with Parameterizations 
Describing Surface OH under Harvard-GEOS Model Conditions 
for July 1994 


Region 

[OHl.xlO 5 

molecules/cm 3 

RMS Error, 

% 

Tropics (30°S-30°N) 

Ocean 

- 10-20 

<8 

Land (desert) 

-20-40 

<8 

Land (forest) 

<10 

8-15 

Land (biomass burning) 

- 10-20 

8-15 

Midlatitudes (30°-60°N) 

Ocean 

- 10-20 

8-10 

Land 

<10 

8-15 

Midlatitudes (30°-60°S) 

Ocean 

<10 

10-30 

Land 

<10 

10-30 

High latitudes (>60°N) 

<10 

10-20 

High latitudes (>60°S) a 

«10 

-- 


a [OHj set to a climatological mean value [Spivakovsky et ai, 1999]. 


A 


V a-x. + a Kr , , > 0 
II /V + 1 

i = 1 

(2) 

N 


y a-x. + a hJ , t < 0 
11 Tv + 1 

1 = 1 

(3) 


Determine coefficients a { (for t = 1, ... , V+l) that define a mini- 
mum of p(a lt ... , a^ + j). We used the simplex method of multidi- 
mensional minimization [Nelder and Mead, 1965; Press et ai , 
1989] to solve for a r Multiple calls to the LSF procedure are nec- 
essary for solving this minimization problem; however, for a sec- 
ond-order polynomial these calls are inexpensive. (To facilitate 
this process, the LSF procedure is interfaced with the point-model 
calculation to generate additional points as needed.) After the 
plane of optimal subdivision is determined, the full procedure 
with selection of appropriate higher-order terms is carried out for 
each subdivision. Construction of polynomials and division of the 
domain are repeated for each subdomain until the accuracy speci- 
fied by the user is achieved, or until the subdomain size is less 
than a user-specified minimum. 

3. Description of Parameterization for OH 

The independent chemical, radiative, and meteorological vari- 
ables chosen to describe parameterized OH are listed in Table 1. 
The choice of independent variables and ranges was guided by the 
current understanding of OH chemistry [e g., Logan et al. , 1981; 
Sillman et ai , 1990; Eisele et al. , 1997; Singh et ai, 1995]. 

Separate parameterizations were obtained for the lower, middle, 
and upper tropospheres (1020-800, 800-300, and 300-150 mbar) 
to account for differences in the chemical regime due to differ- 
ences in the concentration of water, sunlight, and temperature. 
Parameterizations were also obtained for different seasons and lat- 
itude bands because seasonal variations can result in vastly differ- 
ent levels of OH and specific functional dependencies [Kleinman, 
1991]. Due to the varying dependence of OH on the concentration 
of NO,, we constructed separate parameterizations for clean, mod- 
erately polluted, and polluted regions [Lin et ai, 1988]. Finally, 
since variations in isoprene concentrations can act to reduce OH 
concentrations by orders of magnitude [Jacob and Wofsy, 1988], 
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(x 1 O 5 molecu!es / /cm 3 ) (x 1 0 5 molecules/cm 3 ) 

Figure 1 . Parameterized OH (xlO 5 molecules/cm 3 ) versus “true” OH from point-model calculations for four 
locations: (a) clean tropical surface, (b) southern midlatitude surface, (c) tropical surface affected by biomass 
burning, and (d) northern midlatitude upper troposphere. 


we developed parameterizations with and without isoprene. Addi- 
tional domain divisions were determined by the convenience of 
using the resolution of the O 3 and temperature climatologies of 
Logan [1999a,b] and Randel [1992], respectively. In total there 
are about 220 separate parameterizations to describe the tropo- 
spheric OH field. 

Ranges for independent variables used in the parameterization 
are listed in Table 1. Due to a large number of individual subdo- 
mains, it is not feasible to list ranges of independent variables for 
each subdomain. These ranges for chemical species and physical 
parameters were estimated from observations and monthly aver- 
age values predicted by a global model. The global model used 
for estimating ranges was the Harvard-GEOS model of transport 
and chemistry driven by assimilated meteorological data from the 
Goddard Earth Observing System Data Assimilation Office 
(GEOS-DAO) [Bey et ai , 1999]. Since a potential application of 
the parameterization of OH developed here is a study of long-term 
trends of tropospheric CH 4 , we expanded its range to include pre- 
industrial and possible future concentrations. 


4. Results 

RMS errors associated with the parameterizations used to 
describe the surface OH field under Harvard-GEOS model condi- 
tions for a day in July 1994 are shown in Table 2. The RMS error 
given in percent of the mean OH for the sample is used as a mea- 
sure of accuracy. Since removal by OH of CO, CH 4 , and other 
long-lived species occurs predominantly in the tropics [Logan et 
al. t 1981], the highest accuracy was required for 30°S to 30°N. 
We also required higher accuracy in midlatitudes during summer 
because of high concentrations of OH in this season. We tolerated 
lower accuracies where OH concentrations are generally sup- 
pressed due, for instance, to high isoprene emissions. In these 
regions of low OH we chose not to invoke further divisions by the 
parameterization procedure (i.e., to decrease the RMS for the sub- 
domain). In subdomains where OH is very low or negligible due 
to low sunlight (e.g., high latitudes in winter), concentrations of 
OH were set to climatological mean values as a function of lati- 
tude, altitude, and season [Spiv akov sky et ai, 1999]. Lower accu- 
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racies (15-45%) were also tolerated for parameter] zations of the 
upper troposphere (not shown in Table 2) because rates of reaction 
of OH with CO and CH 4 (i.e., two trace gases intended for study 
using the parameterization of OH) are relatively small compared 
to lower tropospheric rates. 

Figure 1 shows parameterized OH versus “true” OH from 
point-model calculations for four subdomains. Since the least 
squares method minimizes the sum of squares of deviations from 
the “true” points, the polynomial fit within a subdomain for higher 
values of OH results in smaller relative errors than for lower val- 
ues. 

5. Summary 

Building on the method of Spivakovsky et al. [1990a], we have 
constructed an up-to-date parameterization of OH which accu- 
rately represents concentrations of OH predicted by a full mecha- 
nism of O 3 , NO fI and hydrocarbon chemistry. This new 
parameterization incorporates current chemical kinetic and photo- 
chemical data, and accounts for reactions of hydrocarbons with 
OH. The computational cost of simulating tropospheric OH is 
reduced by about a factor of 500 when the full chemical mecha- 
nism is replaced by the parameterization of OH. The parameter- 
ization of OH is publicly available from our website 
(http://www-as.harvard.edu/chemistry/trop/index.html). 
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Possible Future Roles 
for 2-D models in GMI 


Malcolm Ko and Charles Jackman 

With inputs from AER, GSFC,a nd LLNL 2-D modeling groups 

Presented at GMI ScienceTeam meeting 
January 4-6 2000 


Possible roles 


• Participate in scientific studies such as M&M, as 
planning tools, as diagnostic tools 

• Serve as platforms for testing components: winds, 
numerical scheme, chemistry, PSC, sulfate 

• Sensitivity of assessment results to different 
transport and chemistry treatments, end to end 

• 2-D transport parameters from 3-D winds. If 
equivalence is established, 2-D model can be used 
to explore parameter space. 




Comparison of Model Transport 
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Grid Resolution 
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Transport Scheme 
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Comparison of Model Chemistry 
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Other issues 


• Long-term averaged (monthly) winds vs. 
“on-line” wind 

• Interannual variabilities 
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Global Implications of the Ozone Loss in a Space Shuttle 

Wake 


M.Y. DANILIN, M.K.W. KO, and D.K. WEISENSTEJN 

AER, Inc., 840 Memorial Drive, Cambridge, MA 02139 


Key Question: 

How important is ozone loss in a Space Shuttle wake on 
the global scale? 


AEAP Annual Meeting, Snowmass, CO, June 2000 



Introduction 

Space Shuttle emits 68 tons of chlorine per launch into the 
stratosphere. 

In situ measurements [Ross et al., 1997;2000] and model 
calculations [e.g., Danilin , 1993; Zittel , 1994] show a severe 
ozone depletion (up to 100%) in a Solid-fueled Rocket Motor 
(SRM) wake if emitted HC1 is converted CI 2 via afterburn- 
ing processes. 

Existing global model calculations [ Prather et ah, 1990; 
Jackman et al . , 1996,1998] of the ozone response to SRM 
launches ignore ozone loss in the SRM wakes. 

Is this assumption justified? 

Addressing this question, we estimate the up- 
per bound of the ozone depletion in the Space 
Shuttle wake and its implications on the global 
scale. To achive this goal, we assume that ALL 
chlorine is emitted as CI2. 



Far- Wake Model 

AER far- wake model [Danilin et al, 1997] with updated 
reaction rates [Sander et al . , 2000] is used. The model as- 
sumes a uniform distribution of species across the wake and 
starts calculations 10 s after exhaust. 

Wake Dilution Scenarios 

R(z, t) = R(t 0 ) + b(t - t 0 )(p(20km)/p(z)) 1/ ' 2 , (1) 

b=1.75 m/s (baseline) and 7 m/s (fast), initial radius R(t o )=200m, 
p(z) is the pressure at altitude z. 



Mixing with Ambient Air 

/ , .r.x n i(z, t)S(z, t) + Ilf (z, t)AS(z, t) 

" i(z - 1 + *> = sSTTJt) • (2) 

n*(z,t) and n?(z,t) are the concentrations of the i-th gas inside 
and outside the wake, S(z,t) is the wake cross section. 




Model Initialization 

The vertical profile of chlorine release from the Shuttle SFtM accord- 
ing to Whitten et al [1975] is used in this study. The main compo- 
nents of the Shuttle SRM are A1 2 0 3 (30% by weight), CO (24%), chlo- 
rine(21%), H 2 (10%), N 2 (9%), C0 2 (4%), and H 2 (2%) emitted up to 42 
km [Prather et al , 1990]. The Table below shows initial concentra- 
tions of Cl 2 in the wake with radius of 200 m. Initial concentrations 
of other exhaust products are obtained from the Cl 2 values by scaling, 
nj (z) =n c j 2 (z) X hr* x fj,d2 / (frcZ2 X /x<) . Other species are taken from the AER 
2-D model climatology at 29°N. 


Vertical profile of CI2, R=200 m 



We split photochemical and wake dilution processes. The ozone re- 
sponses shown are the difference between two runs, with and without 
Shuttle emissions. 



Ozone response to chlorine emissions: HC1 vs CI 2 

Ozone Depletion (in ppmv) for Baseline Dilution, HCI only 



0.01 0.1 1.0 10. 100. 

Time after exhaust, hours 
Ozone Depletion (in ppmv) for Baseline Dilution 


HCI 



0.01 0.1 1.0 10. 100. 

Time after exhaust, hours 



Ozone depletion in the wake is very sensitive to the form 

in which chlorine is released. 



Baseline Case 




Ozone depletion (in % of its background values) and mixing 
ratios of CIO, CI2, and CI2O2 (in ppbv) in the wake for a noon 
Shuttle launch in June. 



Altitude, km Altitude, km Altitude, 


Ozone Sensitivity to Wake Dilution Rate 



Time after exhaust, hours 



Time after exhaust, hours 



Time after exhaust, hours 




From the far-wake to global model 





If [P-L(Oa)] wake and [P-L( 03 )] am fo eni are within ±1%, the 
far-wake model results are ready for use in the AER 2-D 
model. We run the far-wake model 2 days to meet this crite- 
rion. 






Theoretical estimates 

Let consider a box where ozone has a photochemical life- 
time r and is periodically perturbed with amplitude AO 3 
and period T. 



Asymptotic value 

Definition: AO^ = (AO3 + AO^) xe -T / r 
Solution: AO| s = AO3 x 


Definition: 
A Of 


Averaged value 


A Ol + A 01 s 
T 




r 

T 


Steady-state value 

Definition: AO| s = AO™(z, 2 days)S w (z, 2 days)Az 1 ^, 
T=365 days / 12 launches 




Estimates of the global effects of the wake ozone loss 

Knowing ozone loss in the wake, wake size and photochemical lifetime 
of ozone from a global model, one can estimate global effects of the wake 
ozone depletion. 


z, km 

AO™ x VOL, xlO 29 

r, days 

AOf, xlO 29 

10 

0.2 

16.5 

0.11 

12 

1.8 

13.8 

0.82 

14 

2.2 

15.2 

1.10 

16 

2.4 

15.0 

1.18 

18 

10.4 

13.9 

4.75 

20 

12.6 

14.3 

5.92 

22 

14.2 

14.4 

6.72 

24 

15.6 

13.2 

6.77 

26 

18.4 

9.9 

5.99 

28 

•• 26.4 

7.6 

6.60 

30 

44.8 

4.7 

6.92 

32 

62.4 

3.0 

6.15 

34 

87.8 

2.1 

6.06 

36 

* 109.2 

1.6 

5.74 

38 

80.8 

1.3 

3.45 

40 

57.0 

1.0 

1.87 

42 

18.0 

0.8 

0.47 

10-42 

564.2 

- 

70.62 


Peak value: 

5.64 x 10 31 /(305 x 2.687 x 10 26 #/fcm 2 x 3.71 x 10 7 fcm 2 ) = 1.8 x 10“ 3 % 
Steady-state value: 

7.06 x 10 30 / (305 x 2.687 x 10 26 #/&m 2 x 3.71 x 10 7 fan 2 ) = 2.2 x 10“ 4 % 



2-D Model Results 

First run: 12 Shuttle launches per year from Cape Canaveral (29°N) inject 
0.804 kt chlorine into the stratosphere; 

Second run: Additionally, the wake ozone loss is simulated by removing 
AO| I> (z) (in #/km 3 ) between 10 and 42 km altitde in the 23°N-33°N band: 

AOf(z) = AO^(z,2days) SW(z ’ 2 D days) 

^grid 

90°N 

60°N 

30°N 
a> 
xJ 


30°S 

60°S 

90°S 



JFMAMJ JASOND 

Month 



Day of Year 





Ozone Impact: Shuttle vs Aircraft 



Shuttle: Effective Emission Index EEI = -(1-30) 
g(0 3 )/g(fuel) 2 days after exhaust (upper bound!!!). 

HSCT: assuming 1.1 xlO 9 #/cm 3 is depleted in the wake 
with volume of 1 km x20 km 2 due to 8.9 kg of fuel burnt 
[Danilin et al, 1997], one gets negative effective EEI=- 
1 . 1 x 10 24 # /km 3 x 20km 3 x 48 g/ mol/ (6.02 x 10 23 1/ mol x 8900 

g) = -0.2 g(0 3 )/ g(fuel), or an order of magnitude 
smaller than in the Shuttle case at 20 km. 

Subsonic: Ozone production! A0 3 =+0.02 ppbv [Petry et 
al . , 1998], wake volume =50 km 3 , EEI=+0.02 g(0 3 )/kg(fuel) 




SUMMARY 


• Upper bounds of the ozone loss in the Space 
Shuttle wake are estimated assuming all chlorine 
emissions are released as CI 2 ; 

• Local ozone depletion in the Space Shuttle 
wakes could reach 70-80% minutes after exhaust; 

• Global effects of the ozone depletion in the 
wake are larger for faster wake dilution and sum- 
mer sunrise launches of the Space Shuttle; 

• Global effects of the ozone loss in the Shut- 
tle wake are at least an order of magnitude less 
than that due to globally spread chlorine emis- 
sions from the Space Shuttle; 

• EEI(0 3 ) I shuttle ~ 10xEEI(O 3 )|hsct at 20 km (caveats!) 
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